The influence of the weak bond-energy dimerization on the single-particle optical 
conductivity of quasi-one-dimensional systems 
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The single-particle contributions to the optical conductivity of the quasi-one-dimensional systems 
has been reexamined by using the gauge-invariant transverse microscopic approach. The valence 
electrons are described by a model with the weak bond-energy dimerization, while the relaxation 
processes are taken into account phenomenologically. It turns out that the interband conductivity 
of the insulating half-filled case fits well the single-particle optical conductivity measured in various 
CDW systems. In the metallic regime, for the doubled Fermi vector 2fep close to n/a, the conduction 
electrons exhibit a non-Drude low-frequency response, with the total spectral weight shared between 
the intraband and interband channels nearly in equal proportions. For 2fcp — n/a not to small, the 
behaviour of the conduction electrons can be described as the response of a simple Drude metal. 

PACS numbers: 78.20.Bh,71.45.Lr 
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I. INTRODUCTION 

In recent years there has been a lot of experimental 
0, |, |, [§ |, |, § and theoretical § g |Io) interest in 
low dimensional systems which exhibit a gap (for T be- 
low the critical temperature T c ) or pseudogap (T > T c ) 
features in the low-frequency optical conductivity. In 
principle, it is possible to solve the single-particle contri- 
butions to the response functions for the ideal electronic 
system and then to take the single-particle damping phe- 
nomenologically. However, a good agreement with ex- 
periments cannot be achieved by such phenomenological 
extension of the standard transverse microscopic theory, 
which was the question raised by Kim at al. In- 
stead, they developed a more rigorous method which fits 
well the single-particle conductivity in various quasi-one- 
dimensional (Q1D) systems with the CDW instability. 

Surprisingly, a straightforward calculation gives that 
the longitudinal microscopic response theory, with the 
phenomenologically treated disorder, also agrees with 
the single-particle optical conductivity measured in the 
CDW systems The origin of the inconsistency be- 

tween two microscopic approaches represents an interest- 
ing question which will be examined in the present paper. 
The purpose of the paper is threefold; first, to determine 
the gauge-invariance requirement for the electronic sys- 
tems with the bond-energy dimerization (by considering 
the simplest case); second, to argue that in the weak- 
splitting limit the dielectric function has the usual form, 
at variance with the conclusions of the analysis based on 
the electric-dipole approximation |L2| that the dielectric 
function has the Lorentz-Lorenz form even in this limit- 
ing case; finally, to illustrate that in the weak-splitting 
limit various Drude-like and non-Drude low-frequency 
metallic behaviours are possible. 

In this paper, we shall consider the electronic system 
stabilized into an ordered phase described by the wave 
vector Q = n/a. It will be assumed that Q differs from 
the doubled Fermi vector, 2fcp, leading to a metallic state 



with the bond-energy dimerization. Such assumption is 
natural, in particular in the strong-splitting case, describ- 
ing primarily the molecular crystals with an important 
role played by the local field corrections. According to 
Fiefs. |l3|, the weak-splitting Q ^ 2kp microscopic 
model should also be interesting. The optical conduc- 
tivity obtained in this case is expected to be relevant to 
the CDW systems, both at T < T c and T > T c , but for 
the frequencies well above the frequency of the pinned 
collective mode. 

In Sec. 2 we determine the optical conductivity 
a^ t&[ (uj) and the associated dielectric function in the 
Q1D model with the weak bond-energy dimerization, 
taking properly into account gauge invariance of the re- 
sponse theory. Sec. 3 illustrates that the gauge invari- 
ance enables the rigorous treatment of the effective mass 
theorem (and the associated effective number of conduc- 
tion electrons), as well as of the conductivity sum rules. 
Additionally, it is explained how the standard transverse 
treatment of the interband optical processes Jl^, [ll| can 
be remedied. The development of cr* otal (w) with dop- 
ing is shown for a few characteristic cases. Finally, the 
concluding remarks are given in Sec. 4. 



II. THEORETICAL MODEL 

We consider the Q1D electronic system with the bond- 
energy dimerization along the x axis in the presence of 
the electromagnetic fields. The Hamiltonian reads as 



H 
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[tyPia&na + t<a] l+2aa Pna (1) 

+tb(al +brr a n(T + Pl +b Jna) + H.c], 



where the influence of the external fields is introduced 
through 



ft - ;t p ie/(Rc)(R„+ ri )-A(R„+n) 
'ncr l na c - 



(2) 
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and where U na , I S {a, (3}, is a creation operator on the 
R„ + r; site. In order to preserve a simple form of the 
intracell phase factors, we approximate the positions of 
atoms within the unit cell of the dimerized case with their 
positions in the original lattice (for example, r ax = and 
Tpx = a, with ai = ax and a2 = by representing the 
primitive vectors of the original lattice). For the bond 
energies we assume that ti < and |£>| > |i<| > \tb\, 
i.e., the x axis represents the strong axis of the present 
Q1D system. 

In principle, it is possible to generalize the theory to in- 
clude all relevant relaxation processes explicitly. Instead, 
we adopt an approach which assumes that in the clean 
(i.e. underdamped) limit one can separate the relaxation 
of the electron momentum on the static dimerization po- 
tential (or on the (quasi)static Q = 7r/a phonon modes) 
from other electron scattering processes (related to the 
impurities and other phonon modes), calculate the dimer- 
ization (or (quasi)dimerization) effects exactly, and, fi- 
nally, take the remaining relaxation effect into account 
phenomenologically. 



A. Bare Hamiltonian 

The Bloch operators of the dimerized case satisfy the 
relation 

Ll = J= 2^^^. U kx (L,l)lL, 

(3) 

where e %k * Tl * are the bare intracell phase factors, and 
where Uk x (L,l) represent the transformation-matrix el- 
ements dependent on an additional phase tp(k x ), as dis- 
cussed below. For U < 0, the band indices L = S and 
L = A stand for the lower (symmetric) and upper (an- 
tisymmetric) band, respectively. It is also assumed that 
the lower band is partially filled and the upper band is 
empty. If the Fermi level is in the upper band, the fol- 
lowing expressions remain the same, only the electron 
picture has to be replaced by the hole picture. 

Inserting A(R„ + r ; ) =0 in Eqs. (01) and (g), we 
obtain the bare Hamiltonian Hq. After diagonalization, 
this Hamiltonian becomes 

#o = ^E L (k)Lt a L Ua . (4) 

kaL 

The Bloch energies are given by 

£ Aj s(k) = E b (k v )±E a (k x ), (5) 

where 

Eb(k y ) = -2\U\ coskyb, 

E a (k x ) = v/e 2 (M + A 2 (M- (6) 



The corresponding transformation-matrix elements are 
(U kx (S,a) U kx (S,{3) \ 

= ~2~ I e iv(fc * )/2 _ e - i v(fcx)/2 J ' ( 7 ) 

where the auxiliary phase tp(k x ) satisfies the p-symmetry- 
gap relation 

tan^) = (8) 

Other notation is given by 

e(k) = 2\t a \ cosk x a, 
A(k) = A sin k x a, (9) 

with 2\t a \ = \ty \ + |i<| being the averaged intrachain 
bond energy and Ao = \t>\ — \t < | is the magnitude of the 
gap function. 

One recognizes in Eqs. (|^) and (^|) the total intra- 
cell phases of the form & a (k x ) = k x r a + (p(k x )/2, for 
the \an) states, and &p(k x ) = k x rp — ip(k x )/2, for the 
\(3n) states. At this point it is interesting to note that if 
we take the strong-splitting limit, where A /(2|t a |) 1, 
we obtain ip(k x ) ss k x a, and finally <& a (k x ) « ^/3(k x ). 
Such structure of the intracell phases suggests the usage 
of the molecular representation, based on the molecular 
orbitals l/\/2(|cm) ± |/3n)), with a further clcctric-dipole 
approximation (for example, see Ref. 1T3] ) . In the weak- 
splitting limit, however, the atomic |Zn)representation is 
more adequate, as argued in the below discussion of the 
interband current vertices. 



B. Electron-photon coupling 

In order to write explicitly all relevant electron-photon 
coupling functions, we need the Taylor expansion in the 
vector potential of Eq. (|l|) to the second order. The first- 
order term and the direct second-order term become 

kcr LL' i 

x4+ qiI X. + H.c], 

kcr LL' i 

x (10) 

respectively, with X{ S {x,y}. Furthermore, A Xi (q±i) 
and A x .(q±i) represent the Fourier transforms of A Xi (r) 
and A x .(r), and qj_.; ■ sl l = 0. The explicit form of the 
coupling functions, the current vertices J x lL2 (ki) and 

the bare Raman vertices 7 a , 1 a; 2 (fc^; 2), is given in the Ap- 
pendix. 
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It is essential to notice that in the bond-energy dimer- 
ization model two physically different terms comprise 
the interband current vertices. The first term is disper- 
sionlcss and corresponds to the intramolecular interband 
processes (between the states \an) and \j3n)). The sec- 
ond term describes all other, intermolecular interband 
processes. In the strong-splitting limit, one obtains the 
interband current vertices consisting mainly of the in- 
tramolecular contribution, J x A {k x ) m iea\t > \/h, giving 
the associated dipole vertex of the form ea/2. Although 
in this case the response functions are complicated func- 
tions of the short-range dipole-dipole interactions, it is 
possible to apply the electric-dipole approximation, and 
to obtain the Q1D optical model extensively studied by 
Zupanovic at al. |12| . In the present article we want to 
extract physical information from the optical conductiv- 
ity determined in the opposite (weak-splitting) limit. In 
this case the intramolecular and intermolecular interband 
processes are equally important, and, consequently, the 
local field corrections become negligible. 



C. Optical conductivity in the weak-splitting limit 

To obtain the gauge invariant form of the (transverse) 
optical conductivity in the ideal multiband models, it 
is necessary first to gather the diamagnetic photon self- 
energy with the real part of the static interband photon 
self-energies jlTj This results in the effective mass 
theorem [Il0| 



m m 



e 2 E AS {k x ) 



(11) 



Here E AS (k x ) = E A (k) - E s {k) = 2E a {k x ). j^ Xi (k Xi ) 
is the static Raman vertex, which can also be written in 
the form (m/h 2 )d 2 E s (k)/dkl t (T|. When the local field 
corrections are absent for the symmetry reasons, or are 
negligible, the remaining contributions to the total pho- 
ton self-energy lead to the ideal intraband and interband 
optical conductivities given by 



< tra (^i) 
< tor (^ 2 ) 



1 en 



2 - cS 



u> m huj + irji ' 
uV h "ASK 



(12) 



friu - E AS {k Xz ) + ir) 2 
Huj + E AS {k Xi ) + ir]2 



(13) 



vertex over all occupied states [|[ , 

< = ^E^l(k)/5(k), 



(14) 



or as an integral of the square of the electron group ve- 
locity over the Fermi surface, 



,off 



\v x ,(k Xz )\ 2 6[E s (k) - ^ 



(14') 



(notice that for the free electrons both of these expres- 
sions reduce to the total number of conduction electrons) . 
Here /s(k) is the abbreviation for the Fermi-Dirac func- 
tion f\Eg(k)]. Furthermore, notice that the expressions 
( |l2| ) and (Oh describe the case of a normal photon in- 
cidence, A(r) = (A x (y) , 0) , or the case of the parallel 
photon incidence, A(r) = (0, A y {x)). The generalization 
of these expressions is straightforward, and will be briefly 
discussed below. 

It is rather a good approximation to treat the single- 
particle damping (connected with the impurities and the 
remaining phonons) phcnomcnologically, by replacing the 
adiabatic terms rjj with the corresponding damping en- 
ergies £j Xi . Then the total conductivity becomes 



total 



,_intra 



(W, Y,2xi)- 



(15) 



Finally, we obtain that the associated dielectric function 
reads as 



C4 M 



47ri total/, A 



(16) 



The contributions of the high-frequency optical pro- 
cesses, not included in Eq. ([[]), are assumed to be con- 
stant up to frequencies well above the considered fre- 
quency region, and are represented by — 1. 



III. RESULTS AND DISCUSSION 



In the weak-splitting conductivity (pL 5|), 2\t a \ is the 
largest parameter of the model. In spite of the phe- 
nomenological description of the single-particle damp- 
ing, this expression is a rather complicated function of 
the remaining parameters, 2|i&|, k^T, and E A s(kp) 
(the dependence of cr*° tal (o;) on the parameter 2kp — Q 
is represented here through the bare interband threshold 
energy E A s{kp)) which makes further rigorous analytic 
treatment very difficult. In this respect, we shall limit 
the below discussion only on some simple limiting cases. 



A. Strictly ID limit 



respectively. The effective number of conduction elec- 
trons can be shown as an integral of the static Raman 



If the electronic system is treated as strictly one- 
dimensional, i.e. \t{,\ = 0, one obtains that the electron 
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dispersion in the perpendicular direction, Eb(k y ), van- 
ishes. Similarly, for the incidence of the electromagnetic 
fields normal to the chains, the probability of the pho- 
ton absorption (or emission) disappears, as well. In this 
case, the real part of the conductivity comprises only the 
5{oS) term, corresponding to the elastic photon scatter- 
ing. According to the effective mass theorem (fil|), the 
coupling constant arises from two mutually competing 
terms, the diamagnetic term and the static contributions 
in the intcrband term, and is given by the r\\ — ► limit of 
@, with given by (|l4|). For non- normal incidence, 
the effective number of conduction electrons, n| ff , has to 
be replaced by sin 2 ip n° ff , where ip is the angle between 
A(r) and x. In this way, one can derive the well-known 
anisotropic dispersion of the longitudinal plasma modes. 
Namely, 



i e 2 nl s 



sin 2 (p, 



and finally 



^pi(q) 



2^cff 



(17) 



(18) 



To proceed with the numerical calculation, we will take 
the conductivity ( |T5| ) in the usual Q1D limit, \t a \ ^> 
\tf,\ (hereafter this limit is denoted by |t&| — » 0). 
In the optical conductivity considerations concerned only 
with the single-particle contributions, the "explicit" in- 
fluence of the interchain hopping processes on Eq. (15), 
through the Bloch energies, can be neglected. Neverthe- 
less, these processes will be taken into account implicitly, 
since they lead to a finite probability of the photon ab- 
sorption (emission). In such circumstances the ID ver- 
sion of Eq. (|l5|), with the non-zero damping energies, 
can safely be used. Furthermore, we restrict the analysis 
on the case of normal photon incidence, where only the 
Xi = x conductivity is non-zero. We also assume that 
Si = S2 = S and kp = k^x. 



B. Gauge invariance 

The gauge invariance of the conductivity ([D]) can be 
easily proved term by term, by the direct calculation of 
the longitudinal dielectric function, following the proce- 
dure developed in Refs. Q and jl0|]. The most im- 
portant qualitative consequences of the gauge-invariance 
requirement are as follows. 

Fig. 1 shows the influence of the dimerization on the 
(normalized) optical conductivity in the most interesting 
metallic case (/cf = 0.4757r/a), where the TwOK total 
spectral weight is shared between the intraband and in- 
terband contributions nearly in equal proportions. The 
participation of the intraband (i.e. Drude) part in the 
total spectral weight raises significantly with increasing 
temperature. On the other hand, the temperature depen- 
dence of the dc conductivity <r^ c = Re{er* otal (0)} results 
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FIG. 1: The explicit temperature dependence of the real part 
of the optical conductivity for \tt\ — > 0, Ao = 0.5|t a |, E = 
0.2|t a |, and hp — QA75n/a. The temperature scale To is given 
by k B T = 10" 3 • 2\t a \ (for example, T « 4 K, A « 0.1 eV 
and E as 0.04 eV, for \t a \ » 0.2 eV). The intraband and 
interband contributions at T = O.ITq are also shown. 
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FIG. 2: The real part of the optical conductivity as a func- 
tion of the damping energy E, for — > 0, Ao = 0.5|t a |, and 
&f = 0.5ir/a. The filled squares illustrate the prediction of 
the usual, gauge-non-invariant (GNI) transverse model. The 
experimental data measured in K0.3M0O3 (Ref. Q]) are rep- 
resented by the dotted line. 



from two opposite effects, from this explicit temperature 
trend and from the implicit temperature dependence in 
S(T) and/or Aq(T). The experimental data illustrating 
this competition, measured in various low-dimensional 
systems, can be found in Refs. J|, H 0], while a similar 
theoretical analysis is given in Ref. |J. 

More importantly, this figure illustrates that the inter- 
band contribution vanishes for w — > 0, independently of 
kp . Such behaviour is due to the gauge-invariance factor 
[huj / EAs{k x )] 2 in Eq. ([l3]). The physical meaning of this 
observation is that in the metal-to- insulator phase transi- 
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tions (for example, for kp — > 0.57r/a, Ao 7^ 0) the dc con- 
ductivity has to vanish. The usual transverse approach, 
where the gauge-invariance factor is absent, gives, on the 
other hand, a finite crj} c not only in the metallic regime, 
but also in the insulating regime (see Fig. 2 and an exten- 
sive discussion done in Ref. ||). These conclusions are 
valid not only for the bond-energy-dimerization optical 
models, but also for the site-energy-dimerization optical 
models |ll| , as well as for the single-particle contributions 
to <x£ otal in the systems with CDW or SDW instabilities. 

To compare briefly the predictions of the present model 
with experiments, we apply it now to the insulating 
kp = O.bir/a case, and, again, normalize the spectra. For 
the parameter A = 0.5|t o | describing the weak-splitting 
case, we obtain the results, shown in Fig. 2, which fit 
well the single-particle optical conductivity measured in 
the CDW ground state of K0.3M0O3. Notice that the 
usual transverse model (filled squares) significantly over- 
estimates the subgap spectrum. 

Finally, it should be recalled that the gauge invari- 
ant form of CT^, otal allows the consistent treatment of the 
conductivity sum rules ]l0| . For the lower band partially 
filled, the spectral weights of the intraband and total con- 
ductivity are given by 



2 1 'intra 



27re 2 nf. ff 



1 



"total 



m 
2ire 2 1 



kcr 



7 ff(k;2)/ s (k), 



(19) 



respectively, and that of the interband contributions by 
(l/2)(f2 2 otal — f} 2 ntra ). Moreover, the dc conductivity is 
directly related to the intraband spectral weight, accord- 
ing to the well-known relation <j^ c = M} 2 ltra /(47rX;). 



C. \t b \ -> 0, T -> limit 

The characterization of the low-frequency metallic re- 
sponse is intimately connected with the energy region 
between Si and Eas^f) — ^2 (see Fig. 1, for ex- 
ample). In general, in the weak-splitting optical mod- 
els, both the Drude-like behaviour (characterized by 



Re{a x (w)} 



') and the non-Drude behaviour (where 



R-e{CT2,(w)} oc uj~ a ,a < 2) are possible, as can be seen 
from Fig. 3, where the real part of the conductivity for 
various Fermi wave vectors is given. 

A more detailed insight into the low-frequency re- 
sponse can be achieved from the doping dependence of 
n° ff , which is shown in Fig. 4. The following conclu- 
sions are important: (i) The free-electron approxima- 
tion (dashed curve) predicts the electronlikc behaviour 
of the conduction electrons (<9n° ff /<9<5 > 0) in the entire 
doping range < 8 < 2. (ii) The simple tight-binding 
model (dotted curve) leads to the electronlike regime for 

< 5 < 1 and to the holelike regime {dn c f /dS < 0) for 

1 < S < 2. (iii) The dimerized tight- binding model (solid 
curves) is characterized by the electronlike behaviour for 
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= 0.5ll/a,A () = 


— K 


= 0.5*a 


— K 


= 0.475l/a 


k, 


= 0.35ir/a 




1 1.5 

energy / ItJ 

FIG. 3: The real part of the optical conductivity as a function 
of the Fermi wave vector, for Ao — 0.5\t a \ and E = 0.2|£ a |. 
Here crd c (0.57r/a, 0) is the dc conductivity of the &f = 0.57r/a, 
Ao = case. 




FIG. 4: The effective number of conduction electrons Von% 
as a function of the electron doping 6 = 2akp/-K (Vo is the 
primitive cell volume). The prediction of the free-electron 
approximation (FEA) is also given for comparison. 



< 6 < 1 - S c and 1 < 6 < 1 + S c , and by the holclikc 
behaviour for 1 - <5 C < S < 1 and I + S c < 5 < 2. (iv) The 
non-Drude response is expected in the 1 — <5 C <<5<I + < ! > C 
doping region, i.e., where the spectral weight of the inter- 
band contributions becomes comparable to the intraband 
spectral weight. 

For Ao not to large, the critical doping S c is propor- 
tional to Ao. Not surprisingly, for 1 — (5 C < 5 = 1 — 6' < 1 
the electronic system would be described by the effective 
mass approximation which is given by 



-6' = 



in": 



iir; 



(20) 



with 5' and m* x representing, respectively, the number 
of doped holes (with respect to the half filling) and the 
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to the dielectric function in which the local field cor- 
rections are negligible. The analysis reveals two differ- 
ent metallic states. First state corresponds to the usual 
electronlike Drude metal with the optical conductivity 
Re{cr x (u;)} oc w~ 2 , with the effective number of electrons 
characterized by dn x /d8 > and only with one zero- 
point in Re{e^(w)}. Second metallic state is the non- 
Drude state characterized by the Lo~ a behaviour of the 
low-frequency contribution to Re{a x (u>)}, where a < 2. 
There are two zero-points in Ke{e x (uj)}, and the effec- 
tive number of electrons exhibits the holelike behaviour 
dn x s /OS < 0. The comparison with experimental data 
measured in various CDW systems shows good agree- 
ment for the single-particle contributions to the optical 
conductivity. 



FIG. 5: Main figure: The dependence of the real part of the 
dielectric function on the Fermi wave vector, for Ao = 0.5|t a |, 
hy/<LTre 2 /(V m) = 6\t a \, e x = 2.5, and E = 0.2|t„| (see the 
legend of Fig 3). Inset of figure: The influence of the damping 
energies on Ke{e x {^)}, for kp — 0.4757r/a, E = 0.12|t o | (curve 
A), and E = 0.2|t tt | (curve B). 
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effective mass of these holes. It must be noted that the 
effective mass m xx is proportional now to Ao and that 
in the limit Ao — + this mass vanishes. 

Finally, the real part of the dielectric function is af- 
fected by the dimerization in the way shown in Fig. 5. 
In the Drude regime the influence of the dimerization 
is negligible, giving only one zero-crossing of Re{e x (ui)} 
placed at the frequency of the longitudinal intrachain 
(q y = in Eq. (|l8|)) plasma mode. In the non-Drudc 
regime two zeros of the dielectric function appear, de- 
scribing the intrachain plasma oscillations, as well. The 
first frequency corresponds to the intraband plasma oscil- 
lations and the latter one to the total plasma oscillations. 
As illustrated in the inset of figure, in the underdamped 
regime (£, < Ao), considered in this article, the zeros of 
Ke{e x (uj)} are almost independent of the single-particle 
damping. 



APPENDIX A: VERTEX FUNCTIONS 

For the lower band partially filled and the upper band 
empty, relevant are the current vertices (again rp x = a) 



Jy S {k x ) 

T SA 



J x [sin k x a cos (p(k x ) 
A 



2\t a \ 
Jy sin k y b, 



ik x awx<p(k x )], 



Jx (ky) = iJ x [sin k x a sin tp(k x ) 



Ao 
~2|t„l 



cos k x a cos <p(k x )], 



Jy (k y ) 



0. 



(Al) 



(A2) 



and the bare Raman vertices 



IV. CONCLUSION 

In this paper we have studied the influence of the weak 
bond-energy dimerization on the optical properties of the 
Q1D systems, for the electron doping < S < 1. It is 
found that the gap parameter has the p x symmetry. It 
is shown that the weak bond-energy dimerization leads 



7xx( fc *;2) = l xx [cos k x a cos (p(k x ) 
, Ao 



2 It-, 



■ sinfc x asin</2(fc x )], 



7™(fc y ;2) = y° cos kyb. 



(A3) 



(A4) 



Here J° s = eh/ (aim XzXi ), j° iX . = m/m XzXt and m XiXi 
h 2 /{2\t l \a 2 l ) QU\ € {\ta\,\h\}, a, £ {a,b}). 
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